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Abstract 


The  effect  of  electromagnetic  reaction  to  relaxation  of  a  molecule 
adsorbed  on  a  rough  surface  is  considered.  We  include  the  reaction  fields  due 
to  the  system  dipole  relaxing  by  both  photon  radiation  and  decay  to  delocalized 
surface  plasmons  in  the  determination  of  absorption  by  the  adsorbate.  Surface 
plasmon  damping  is  shown  to  be  important  for  both  small-  and  large-sized 
roughness  features . 


PACSNos.:  42 . 50kb ,  68.35.Bs,  82. 65. Pa 


1 .  Introduction 


Following  the  discovery  of  surface -enhanced  Raman  scattering  (SERS) , 

there  has  been  a  surge  in  studies  of  surface- induced  linear  and  nonlinear 
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optical  phenomena.  Among  them,  surface- induced  photochemistry  has  received 

a  great  deal  of  attention.  In  a  photochemical  process,  the  primary  step  is 

absorption  by  molecules  from  an  incident  laser  field  and  consequent  excitation 

to  a  higher  electronic  state,  followed  by  either  dissociation  or  relaxation  to 

a  lower  state,  which  may  include  both  radiative  and  nonradiative  mechanisms. 

The  relative  rate  of  dissociation  determines  whether  the  molecules  undergo  any 

significant  chemical  decomposition.  As  in  SERS,  the  presence  of  a  nearby 

surface  modifies  the  local  electric  field  and  thus  affects  both  the  primary 

absorption  and  the  relaxation  steps. ^  ^  The  increase  in  the  local  electric 

field  is  understood  to  be  due  to  resonant  excitation  of  radiative  plasma 
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oscillations  of  the  surface.  ’  ’  ’  There  may  be  additional  competing,  but 
mitigating,  line -broadening  effects  due  to  energy  transfer  to  the  nearby 
surface  from  the  excited  molecules.^ 

Electromagnetic  calculations,  in  the  context  of  SERS,  for  a  single 
molecule  near  a  small  metallic  particle  or  a  rough  surface  show  chat  the 
transfer  of  energy  between  the  molecule  and  the  surface  takes  place  via 

7  8 

excitation  of  electronic  modes  of  the  associated  particulate  structure. 
Radiative  decay  of  these  electronic  excitations  restores  partially  the  energy 
lost  by  the  molecule.  However,  in  the  case  of  a  rough  surface,  the  remaining 
energy  is  gone  permanently  by  resistive  losses  in  the  particle  and  by  the 
resonant  coupling  of  the  electronic  modes  to  the  delocalized  surface  plasmon 

9 

modes  of  the  underlying  substrate.  In  the  case  of  a  single  ellipsoidal  bump 

9 

protruding  out  of  a  semi- inf inite  flat  plane  as  the  model  for  a  rough  surface, 
the  latter  coupling  mechanism  and  corresponding  energy  loss  have  relatively 
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large  races  (an  order  or  two  larger  in  magnitude  compared  to  the  radiative  and 

resistive  decay  rates)  for  large-sized  bumps. 

In  view  of  this,  we  have  studied  the  relative  effects  of  including  the 

reaction  fields  due  to  the  decay  of  the  molecule -bump  system  to  surface 

plasmons  as  well  as  radiation  damping  in  the  electromagnetic  calculation  of 

absorption  near  a  rough  metal  surface.  The  importance  of  the  reaction  field 

due  to  photon  radiation  by  a  system  consisting  of  a  single  molecule  near  an 

12 

ellipsoidal  particle  was  first  recognized  by  Wokaun  et  al  in  the  context  of 
SERS.  The  primary  effect  of  this  radiation  damping  is  a  decrease  in  the  local  - 
field  enhancement  and  hence  a  decrease  in  the  absorption  enhancement  for  large - 
sized  particles.  Because  the  decay  to  delocalized  surface  plasmons  is  the 
dominant  mechanism  for  large-sized  bumps , ^  a  molecule  in  the  vicinity  of  the 
bump  would  experience  the  effect  of  this  decay  channel  via  the  reaction  field. 
It  is  the  intent  of  this  paper  to  study  the  effects  of  these  two  reaction 
fields,  corresponding  to  photon  radidation  and  decay  to  surface  plasmons  on 
absorption  enhancement. 

2 .  Model  System 

A  popular  model  for  a  rough  surface  is  a  hemispheroidal  bump  protruding 

out  of  a  semi- infinite  flat  plane.  The  system,  consisting  of  a  single 

polarizable  molecule  located  at  a  distance  H  from  the  bump  surface  and 

irradiated  by  an  incident  laser  field  Eg,  has  been  treated  in  the 

literature. In  spheroidal  coordinates  (£ ,  7,  ^) ,  the  hemispheroid  with 

aspect  ratio  a/b  is  characterized''  by  a  shape  parameter  £g  -  a/f,  where  f  -  (a 
2  4 

-  b  )  The  material  characteristics  of  the  substrate  are  accounted  for 
through  its  dielectric  function  «(w). 


— t  -4  A  , 

The  incident  laser  field  Eq  is  taken  along  the  z-axis,  Eq  -  EqZ ,  so  that 

the  induced  dipole  m  and  local  field  E^oc  are  aligned  along  the  z- direction. 
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E^qc  consists  of  several  parts  and  is  given  by  ’ 

hoc  -  E0  •  <Vf>  }  Bn  %  «1>  ■  <U 

n 

where  -  (a+H)/f,  are  the  Legendre  functions  of  the  second  kind,  and  the 
"prime"  indicates  the  first  derivative  with  respect  to  the  argument  The 

incident  field  Eq(w)  induces  a  dipole  moment  in  the  molecule  and  also  polarizes 
the  bump-plane  system,  which  in  turn  produces  an  electric  field  E^  at  the 
location  of  the  dipole.  The  near  field  of  the  dipole  polarizes  the  bump-plane 
system,  creating  an  "image"  dipole  which  generates  a  field  E^  at  the  molecular 
dipole  location.  ’  Both  and  Ej  are  contained  in  the  second  term  on  the 
right-hand  side  of  Eq.  (1). 

The  molecule-bump  system  may  be  imagined  as  a  single  dipole,  of  moment  D, 

2 

located  near  a  flat  surface,  where  D  -  fi  +  B^f  /3.  The  excited  system  dipole  D 

then  relaxes  radiatively  by  emitting  photons  and  nonradiatively  by  coupling  to 

the  delocalized  surface  plasmons  of  the  flat  substrate.  The  power  dissipation 

into  either  of  the  above  decay  channels  have  been  calculated  by  Das  and 
9 

Gersten. 

We  account  for  the  energy  loss  of  the  system  by  adding  the  radiation 

reaction  fields  E  .  (for  photon  emission)  and  E  (for  coupling  to  surface 
raa  sp 

plasmons)  to  the  local  electric  field.  To  maintain  self  consistency  we  write 


Eioc  “  Eo  *  fl-a 


(2) 


and 


B  -  B(„,E')  , 


where  -  EQ  +  Erad  +  E  ,  qn  -  (1/f  ,  and  I  represents  a  formal 

solution  for  the  coefficients  B  in  vector  form. 

n 

The  reaction  fields  associated  with  both  photon  radiation  and  delocalized 
surface-plasmon  excitation  by  the  system  dipole,  D  -  D(/*,Eq),  may  be  obtained 
from  the  power  radiated  to  photons  (P  )  and  power  lost  to  the  excitation  of 
surface  plasmons  (P  )  by  the  system  dipole.  From  Ref.  9, 


r/J  JtWtu 

ra  c  1  J0  (e-1)  [  ( 


e+1)  cos  5-1] 


-]  ,  (4) 


p  -  w  |d|2(«)3  - *f/2 - 

SP  C  (-e-iyU(l-e) 


where  u,  the  frequency  of  incident  photons,  is  taken  to  be  resonant  with  the 

surface-plasmon  excitation  frequency.  The  fields  Erad  and  E  are  determined 

by  using  the  expressions  for  power  loss  given  in  Eqs.  (4)  and  (5)  and  demanding 

that  the  work  done  by  the  radiation  reaction  force  on  the  molecular  dipole  is 

12 

equal  to  the  negative  of  the  energy  radiated  during  a  given  time.  This 
calculation  yields 


(e-l)[(e+l)cos  5  -1] 


0 


and 


E 

sp 


4iriD(-)3 

c 


-Lii ). 


(-€-1) 


5/2 


3 

<1-0 


(7) 


For  a  Drude-type  molecule  driven  by  the  local  field  E^oc>  the  power 
absorbed  is  given  by 


P  -  -  “  Im  (n*  Kx  Eq) 


(8) 


where 


2 

w0a0^1^0 

M  -  —  2  (9) 

(wQ-w  -ia>7) 

is  a  self-consistent  expression  for  the  induced  dipole  moment  of  the  molecule. 
Part  of  the  local  field  goes  to  shift  and  broaden  the  excited  level  of  the 
molecule  and  appears  through  the  modified  width  y  -  +  7S •  In  Eq.  (9),  Wq, 

and  7q  are  the  natural  frequency,  static  polarizability  and  natural  line 
width  of  the  molecule  in  the  absence  of  the  surface,  respectively.  The  power 
Pq  absorbed  in  the  absence  of  the  surface  is  obtained  by  taking  -  1  and  y  - 
7q  in  Eqs.  (8)  and  (9).  The  ratio  RA  -  P/Pq  will  be  called  the  enhancement 
factor  for  absorption  near  a  rough  surface . 

3.  Results 

Figure  1  shows  the  absorption  enhancement  ratio  for  a  molecule  located 
40  a.u.  from  the  surface  of  a  silver  hemispheroidal  bump  on  a  flat  silver 


surface  (the  model  rough  surface)  as  a  function  of  the  semimajor  axis  a  of  the 
hemispheroid  for  various  aspect  ratios  a/b.  Optical  constants  for  silver  were 
obtained  from  Ref.  13.  Numerical  calculations  are  done  for  the  molecule 
resonating  with  both  the  incident  laser  frequency  w  as  well  as  with  the  grounch- 
state  resonance  frequency  of  the  bump.  For  a/b  -  2,3  and  4,  the  ground  states 

9 

correspond  to  Ree(u)  -  -7.34,  -11.3  and  -15.9,  respectively.  The  frequencies 

associated  with  the  ground  states  are  obtained  by  fitting  the  optical  data  of 

2 

Ref.  13  to  a  free-electron-like  expression  Ree(w)  -  A  +  B/w  .  Other  molecular 

.  3 

parameters  used  in  the  numerical  analysis  are  y^  -  10  a.u.  and  -  10  a.u. 

These  graphs  show  the  effects  of  both  including  and  excluding  the  reaction 

fields  in  the  calculations.  Let  us  denote  the  contributions  to  the  reaction 

fields  due  to  photon  radiation  and  decay  to  surface  plasmons  by  the  system 

dipole  D  by  urad  and  u  ,  respectively.  The  dashed  curves  are  obtained  for  the 

case  where  the  reaction  contributions  are  neglected,  i.e.,  u  -  urad  +  u  -  0, 

and  the  solid  curves  are  for  u  *  0,  For  a  given  bump,  (a/b  fixed)  the 

absorption  enhancement  with  u  i*  0  is  substantial  for  small  values  of  a,  but 

rapidly  decreases  to  values  below  those  for  u  -  0  with  increasing  a,  with  a 

critical  a  at  which  both  become  equal.  This  is  understandable  since  the 

effects  of  including  the  reaction  fields  (u  **  0)  are  significant  for  large- 

12 

sized  roughness  features  (a  large  with  a/b  fixed) . 

As  mentioned  earlier,  u  consists  of  two  parts:  urad  corresponds  to  the 

inclusion  of  photon  radiation  reaction  and  u  to  the  inclusion  of  reaction  of 

sp 

the  decay  of  the  system  dipole  to  delocalized  surface  plasmons,  which  is  the 

primary  relaxation  mechanism  very  close  to  the  surface.  The  effect  of  only 

u  ,  *  0  was  discussed  in  Ref.  12  in  the  context  of  SERS.  Curve  c  in  Fig.  2 
rad 

represents  the  effect  of  including  only  urad  in  the  determination  of  absorption 
enhancement.  This  is  in  qualitative  agreement  with  the  results  of  Ref.  12. 
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For  smaller-sized  bumps,  the  effect  of  adding  radiation  damping  is  not 

significant,  whereas  for  larger-sized  bumps  the  enhancement  is  limited  by  the 

addition  of  the  radiation  damping  field  in  the  local  field.  Since  we  are  not 

comparing  our  results  to  any  realistic  experimental  observation,  we  cannot  draw 

a  quantitative  conclusion  except  to  say  that  the  qualitative  trends  are  in 

12 

agreement  with  previous  results  in  the  literature.  On  the  other  hand,  for 
optical  absorption  by  a  molecule  near  a  rough  silver  surface,  the  effect  of 
including  the  surface  plasmon  reaction  field  is  substantial  for  both 
intermediate  and  larger- sized  bumps.  This  is  shown  in  Fig.  2  (curves  a  and 
c) .  It  clearly  illustrates  the  importance  of  taking  into  consideration  the 
reaction  field  corresponding  to  excitation  of  surface  plasmons,  which  is  the 
main  thrust  of  this  paper.  It  is  not  our  intention  to  project  curve  a  in  Fig. 

2  to  be  monotonically  increasing  with  decreasing  a.  A  portion  of  it  (for 
values  of  a  <  300  a.u.)  is  simply  not  shown  for  the  relative  scales  involved  in 
order  to  provide  a  qualitative  comparison  among  the  three  curves.  In  fact,  the 
full  nature  of  curve  a  (Fig.  2)  is  depicted  in  Fig.  3  (curve  d) ,  which  shows 
the  quenching  effect  of  ugp  **  0 . 

In  Fig.  3,  we  plot  as  a  function  of  a  for  various  fractions  of  u 

included  in  u.  We  write  u  -  u  ,  +  xu  where  x-0.2,  0.3,  0.5  and  1.0 

rad  sp 

corresonding  to  curves  a,  b,  c  and  d,  respectively.  Aside  from  the  quenching 
in  absorption  observed  with  increasing  x,  there  is  a  movement  of  the  peak 
position  to  smaller  a-values.  This  gives  a  better  understanding  of  the 
enhancement  of  optical  absorption  near  a  rough  surface.  If  one  wants  to 
produce  a  predetermined  absorption  enhancement  (consequently  a  similar 
photodissociation  enhancement)  in  a  molecule  held  at  a  fixed  distance  from  the 
surface,  the  surface  morphology  has  to  be  designed  accordingly.  In  other 


words,  the  surface  preparation  is  of  paramount  importance.  On  the  other  hand, 
if  the  surface  roughness  is  specified,  one  needs  to  put  the  molecule  at  a 
certain  distance  away  to  obtain  the  maximum  enhancement.  For  accurate 
prediction  of  spectroscopic  properties  of  adsorbates,  one  has  to  account  for 
the  reaction  field  due  to  surface  plasmon  excitation  within  the  theory.  Given 
the  surface  preparation  techniques  known  to  date,  it  is  perhaps  possible  to 
test  the  predictions  made  here  regarding  inclusion  of  reaction  fields  due  to 
surface-plasmon  excitation  in  the  calculations. 

4.  Conclusion 

In  conclusion,  we  have  considered  a  correction  to  the  electromagnetic 

theory  of  surface -enhanced  absorption  by  taking  into  account  the  reaction 

fields  due  to  both  radiation  and  decay  to  delocalized  surface  plasmon.  We 

chose  a  model  where  decay  to  delocalized  plasmons  of  the  flat  surface,  besides 

photon  emission,  by  the  system  is  possible.  It  is  found  that  for  small-sized 

roughness  features,  the  radiation  damping  correction  is  unimportant,  but  the 

surface  plasmon  reaction  is  important  for  both  small-  and  large -sized  features. 

It  should  be  pointed  out  that  the  rough  surface  considered  is  only  a  model 

surface.  Realistic  rough  surfaces  are  quite  different  and  difficult  to  model. 

Predictions  of  this  calculation  may  be  checked  experimentally  on 

14 

microlithographically-prepared  surfaces . 
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u  r*  0  and  the  dashed  curves  for  u  -  0.  The  resonance  frequencies 
for  the  ground- state  resonances  of  the  silver  bumps  of  aspect  ratios 
a/b  -  2,  3  and  4  are  0.1,  0.086  and  0.076  a.u. ,  respectively.  The 
molecule -spheroid  distance  is  H  -  40  a.u. 

Figure  2.  RA  as  a  function  of  a  for  fixed  a/b  -  3.  Curves  a  and  b  correspond 

to  u  r*  0  and  u  -  0,  respectively.  Curve  c  is  obtained  when  only  the 

photon  radiation  reaction  field  is  taken  into  account  (u  ,  0,  u 

r  Eu  sp 

-  0).  The  molecule -spheroid  separation  is  H  -  40  a.u.,  e(wQ)  - 
-(11.3)  +  i(0 . 24)  and  wQ  -  0.086  a.u. 

Figure  3.  R^  as  a  function  of  a  for  fixed  a/b  -  3.  Curves  a-d  correspond  to  x 
-0.2,  0.3,  0.5  and  1.0,  respectively,  where  u  -  urad  +  xusp.  The 
other  parameters  are  the  same  as  in  Fig.  2. 
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